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We present a microscopic study of the nuclear magnetic relaxation rate 1/Ti based on the five- 
orbital model for iron-based superconductors. We mainly discuss the effect of the "inelastic" quasi- 
particle damping rate 7 due to many-body interaction on the size of the coherence peak, for both 
s ++ and s±-wave superconducting states. We focus on Ba(Fei_ ;c Co :E )2As2, and systematically 
evaluate 7 in the normal state from the experimental resistivity, from optimally to over-doped 
compounds. Next, 7 in the superconducting state is calculated microscopically based on the second 
order perturbation theory. In optimally doped compounds (T c ~ 30 K), it is revealed that the 
coherence peak on \/T\T is completely suppressed due to large 7 for both s++ and s±-wave states. 
On the other hand, in heavily over doped compounds with T c < 10 K, the coherence peak could 
appear for both pairing states, since 7 at T c is quickly suppressed in proportion to T c 2 . By making 
careful comparison between theoretical and experimental results, we conclude that it is difficult to 
discriminate between s++ and s±-wave states from the present experimental results. 



I. INTRODUCTION 

Since the discovery of superconductivity in 
LaFeAsOi-zF^ji the superconducting mechanism 
and paring symmetry had been discussed intensively. 
In many optimally-doped compounds, the super- 
conducting gap is fully gapped as reported by the 
penetration depth measurement 2 and the angle-resolved 
photoemission spectroscopy (ARPES)i^ As for the 
superconducting mechanism, the s-wave pairing with 
sign change of the order parameter between the hole 
and electron Fermi pockets, so called s±-wave, medi- 
ated by the anti-ferromagnetic fluctuation had been 
proposed from the early stage as a possible pairing 
state in the iron pnictides£~— It is supported by quasi- 
particle interference analysis by Scanning Tunneling 
Microscopy/Spectroscopy (STM/STS) measurement by 
Hanaguri et al. 10 However, the small T c -suppression 
against nonmagnetic impurities is not consistent with 
the s±-wave stated 1 - - — 

On the other hand, orbital fluctuation mediated s- 
wave superconducting state without sign reversal 
wave state) had been investigated. 16 We have shown 
that strong ferro- and antiferro-orbital fluctuations de- 
velop due to the combination of Coulomb and e-ph 
interactions *2r— Consistently, large softening of the 
shear modulus Cqq£-~— and renormalization of phonon 
velocity^ are observed well above the orthorhombic 
structure transition temperature T s . These phenom- 
ena strongly suggest the existence of strong ferro-orbital 
(charge quadrupole O x 2_ y 2) fluctuations, considering the 
large strain-quadrupole coupling.— ~— In addition, exper- 
imental "resonance-like" hump structure in the neutron 
inelastic scattering is well reproduced in terms of the 
s ++ -wave state. 26 i 27 

The nuclear-magnetic resonance (NMR) and the 
nuclear-quadrupole resonance (NQR) measurements are 



useful for the discussion of the pairing symmetry. The 
coherence effect of superconductivity, appearing as a 
Hebel-Slichter peak (coherence peak) in the nuclear spin 
relaxation rate (1/Ti), was one of the crucial experi- 
mental proofs of the Bardeen-Cooper-Schrieffer (BCS) 
theory, characterized by conventional s-wave Cooper 
pairs with an isotropic gapj2£ In the Fe based supercon- 
ductors, many experimental results of 1/Ti have been 
published^ - — It was found that the coherence peak in 
1/Ti is absent in many compounds, like electron-doped 
Ba(Fei_ :E Co a ;)2As2 and hole-doped Bai_ a; K :c Fe2As2. 

The size of the coherence peak had attracted great at- 
tention to distinguish between s++ and s±-wave states: 
In a simple BCS theory, the peak size is larger in the 
s ++ -wave state, while it is reduced in the s±-wave state. 
However, the coherence peak is suppressed by the "inelas- 
tic" quasi-particle damping rate 7, which is prominent in 
moderately and strongly correlated systemsj 36 i 37 For this 
reason, the coherence peak is not observed even in several 
conventional s-wave superconductors with T c > 15 K, 
such as boron carbide YM2B2C 38 and A-15 compounds 
V3SL 39 Thus, inelastic scattering due to many-body ef- 
fect has to be taken into account for a quantitative anal- 
ysis of 1/Ti. 

Impurity effect (— "elastic" quasi-particle damping ef- 
fect) on T c and 1/Ti is also important in studying the 
pairing symmetry. In the s++-wave state, both T c and 
1/Ti are insensitive to impurities. In contrast, the s±- 
wave state is easily suppressed by impurities like the d- 
wave state, according to the study based on the five- 
orbital models The impurity-induced gapless state in 
the s±-wave state would give strong influence on 1/Ti 
for T <C Tc^ i 40 ' 41 However, impurity- induced gapless 
state is realized only when T c is strongly suppressed to be 
~ T c o/3 theoretically. 42 Thus, gapless behavior of 1/Ti 
for T <s^ T c observed in some compounds cannot be ex- 
plained by this scenario. 

Recently, authors in Refs. and HH had shown 
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that the coherence peak in the s++-wave state disap- 
pears when inelastic scattering 7 at T c is as large as 
the superconducting gap at T = 0. However, quanti- 
tative estimations of 7 and its T and u dependences are 
still lacking. Interestingly, recent NMR measurement re- 
ports a small coherence peak of 1/T\ in the heavily over 
doped LaFeAsOi-zF-r with low transition temperature 
T c ~ 5 Ki2£ Now, microscopic study of 1/Ti by including 
inelastic scattering effect is desired to discuss the super- 
conducting pairing state as s± or s ++ . 

In this paper, we investigate the effect of inelastic scat- 
tering rate 7 on the nuclear magnetic relaxation rate 1 jT\ 
for both s++ and s±-wave states. For a quantitative dis- 
cussion on the existence of the Hebel-Slichter coherence 
peak, we employ the two-dimensional five orbital model. 5 
Here, 7 is the key parameter of the present study: At 
T c , the value of 7 is carefully estimated from the ex- 
perimental conductivity^ and its temperature and the 
frequency dependences below T c is obtained by the sec- 
ond order perturbation theory. All results shown below 
are obtained with band filling n — 6.1. 

The contents of this paper are as follows; In Sec. HU we 
explain how to calculate the nuclear magnetic relaxation 
rate \fT\T and quasi-particle damping 7. The obtained 
numerical results are explained in Sec. IIIII Finally, we 
make comparison between theoretical and experimental 
results in Sec. IIV1 and discuss the possible pairing sym- 
metry. 



II. FORMULATION 



A. Green function 



Now, we study the 10 x 10 Nambu BCS Hamiltonian 
%k composed of the five orbital tight binding model 5 - and 
the band-diagonal SC gapjA 5 - The Hamiltonian is given by 



£ = £4^ fe , 



where c\, and dk are vectors, 



at _ (J 



(1) 



(2) 



and cjj k (c at k,a) is a creation (annihilation) operator 
of an electron for band a with wave vector k and spin a. 
Ti.^ is 10 x 10 Nambu BCS Hamiltonian given as 



T/0 



Hk Afc 
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where Hk and A*, are 5x5 matrices in the band-diagonal 
basis, 
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FIG. 1. Fermi surfaces of the 5 bands model for iron pnictide 
in the unfolded Brillouin zone. 



where e Q .fc and A Qj fe are a dispersion and a gap func- 
tion of quasi particle for a and k, respectively. Figure Q] 
shows the Fermi surface (1 Fe atom/unit cell) for electron 
doped iron pnictides. Fermi surface ct\ and around T 
point are hole pockets and a% and 014 are electron pock- 
ets. Hereafter, we approximate A Qi fc as follows: 



A, 



A r 



A° tanh 




(5) 



where A° is a superconducting gap at zero temperature 
and, in the preset work, we assume that the supercon- 
ducting gaps for each band are fc-independent. 

The 10 x 10 Green's function Qk in the Nambu repre- 
sentation is given by^ 5 - 

-1 



Gk{iu n ) = (i 



a - m 



(6) 



where u n = irT(2n + 1) is the fcrmion Matsubara fre- 
quency. The normal self energy Sfe(zw n ) represents the 
inelastic quasi-particle damping and mass enhancement 
due to many-body interactions. We neglect the impurity 
induced self energy since it does not change the density of 
states (DOS) in a simple s-wave state, Eq. (J5J) (Anderson 
theorem). When Efc(iw„) has band diagonal form, Eq. 
© also becomes a band diagonal form, and the Green's 
function for band a is given by a 2 x 2 matrix, 



1 



\Ga,k( iuJ n) 



£a,k + ^a,k(~ 

A„ 



■iu) n ) 



-A a 



ioj n — e a ,k — S Q .fc(iw„) 



G a ,k(iw n ) 
F a ,k{iu n ) 



F a> k(iuj n ) 

-Ga,k{ — iw r , 



(7) 
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where \G a ,k{iw n ) | is the determinant of the inverse ma- 
trix of the Green's function, and G a . k , F ak and E aj fc 
are the normal Green's function, the anomalous Green's 
function and the normal self-energy for band a with wave 
vector k, respectively. The local Green's function g a is 
obtained by the k summation; 



g a (iuj n ) = — ^ Q a ,k(iUn) 



deN(e)G a (e, iu n ), 



(8) 



where N(e) is the quasi-particle DOS. In this e integral, 
the dominant contribution comes from the explicit e term 

of the denominator |5 Q ,fc(«w„) |. By neglecting the e de- 
pendence of N(e) and applying the infinite dimensional 
approximation for E a (e, ioJ n ) (see the following subsec- 
tion ITTCj, we obtain 



-7rAT(0) / iw a (iw n ) A a 
\Ju) a (iu n ) + A z a K 

g a {iw n ) /«(«*>„) 



(9) 



From the five orbital model, the total DOS per spin is 
iV(0) = 0.66 eV" 1 at Fermi level. Then, analytic contin- 
uation (iuj n — > u> = u> + iS) yields the retarded Green's 
function as follows, 



3a M = 
/ Q R (^) = 



-7riV(0)6£M 
-nN(0)A a 



(10) 



Al 



where u) r (uj) is defined as 

w£(w) = u - iImE R (u;) 

= ^ + *7aM, 



(11) 



where 7* is the "renormalized" quasi-particle damping, 
which is described by using "unrenormalized" quasi- 
particle damping 7 = — ImE R and bare and effective 



masses m and m* as 



The mass enhancement factor m*/m is given by 



m 
m 



lim 



<9ReE R (u;) 



(12) 



(13) 



which is reported as 1 ~ 3 for iron pnictides called 
1111 and 122 systems from various experiments such as 
de Haas- van Alphen measurements optical spectral 
weight,— and Seebeck effect with specific heat^S The 
mass enhancement factor has a relatively large value in 
the optimally doped systems and gradually decreases by 
carrier doping. 



B. Nuclear Magnetic Relaxation Rate 

The nuclear magnetic relaxation rate \/T\T in the su- 
perconducting state is given by the standard formula: 



— — cx — lim Im 



(14) 



where Xq ( w ) is the superconducting spin susceptibility 
and given by 



1 I { , UJ + Oj' Uj' 

dui tanh — tanh — — 

V 2T 2T 



x (lmGl k+q (uj + u/)ImG R fe (u/) 

+ImF* k+q (u + c/)ImF£ fe (u/)) . (15) 

In Eq. (|14p , the limitation becomes a differential of hy- 
perbolic tangent and each summation of a, j3, k and q 
can be calculated independently. Then, 1/T\T becomes 

1 [°° , Im 5 R (w) 2 + Im/ R (w) 2 

oc / aw = ^ 

TiT J_ x 



47rTcosh^(w/2T) 
du)X(uj), 



(16) 



with 



(17) 



a 

where the integrand of uj was defined as X(oj). 

C. Second order perturbation theory 



A quasi-particle excitation is well defined if the damp- 
ing rate 7 is small compared to the energy scale. In 
this case, the quasi-particle damping rate 7 can gener- 
ally be computed from the imaginary part of the normal 
self energy^ 



ImE R (a;) = -7fc(w). 



(18) 



By using the second order perturbation theory, we de- 
rive the microscopic picture of 7 and analyze the NMR 
experimental result. We calculate the electron-electron 
scattering based on the single band Hubbard model for 
simplicity. From the second order perturbation theory, 
the self energy E*. due to the electron-electron scattering 
is given by 



T/2 rp 

Efc(tw„)- "'' ^ 



N 



k' ,n' 



where V c g is the effective electron-electron interaction en- 
hanced by spin and orbital fluctuations. Its value can 
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be estimated from the observed conductivity in subsec- 
tion QTD] The imaginary part of the retarded self energy, 
which is obtained by the analytic continuation, is given 
by, 



V 2 
2nN 



E 



doj' [ tanh — — 
, V 2T 

o,r 



coth ■ 



2T 



(20) 



where the superscript "0" on the Green function and sus- 
ceptibility means the absence of self-energy correction. 

In the infinite dimensional approximation, the each 
summation of wave vectors can be taken independently 
and we obtain the quasi-particle damping rate due to the 

I 



electron-electron scattering as 

' - N - ,.i i = - 

N 



7(w) 



N E 

k 



k 



V 2 



— U 



tanh — — 
2T 



coth • 



2T 



xIm 5 > R (w')Imx°' R (w 
with bare susceptibility, 
1 



(21) 



Imx u ' K (w) 



2tt 



, I OJ + OJ CJ 

duj tanh — tanh 



— oo \ 
0,R 



(lm.g u ' H (w 
+Im/°' R (w 



2T 

+ uj')Imf> K (oj')) 



2T 



(22) 



To see the role of coherence factor, we rewrite the imaginary part of the Green's functions with ImG^' R (w) 



-f{(l" 



j*-)6(oj-E k ) + (l 



-£ fc )} and ImF^ R (w) 



_ 7rA fc 
2 -Eft 



{c5(w-£ fc )-5(w + £; fc )}, where £ fe = 



By neglecting the energy- dependence of the DOS near the Fermi level, we obtain another notation of the self energy 
as, 



Im£ R H 



2T 



V 2 cosh ■ 

fe fe cosh cosh cosh 

2T 2T 2T 



1 



A,> ' 1 Aa - ) {*(|w| - |£fe< + Sfc-fc'+fc- - E k „\) + 6{\oj\ - \E k , - E k _ k , +k „ + E k „\)} 

{S(\oj\ - \E h ' - Eh-k'+k" - E h "\) + S(\u\ - \E k , +E k - k , +k „ + E k „\)} [23) 



E k -k'+k"E k 
A k - k > +k » Afc 



E k - k >+ k "E k n 



For T <T C , the first term with the coherence factor (1 + AA' /EE') gives the dominant contribution for the s_|__|_-"W"ave 
state. At T = 0, because of the thermal factor, the self energy becomes 



8A 2 



Ek-k'+k"E k r. 



S(\oj\ 



\E k > + -Efc-fc'+fc" + Sfe// 



(24) 



It is clear from the delta function that it equals for u < 3A m ; n . 

I 



Experimental studies on the iron pnictides such as 
ARPES 52 , point contact tunneling, 53 NQR 31 and spe- 
cific heat 5 -^ have demonstrated that there are at least two 
different superconducting gaps, small gap As and large 
gap Al. Multi-gap superconductivity has been seen in 
a number of systems including MgB 2 which is an s-wave 
superconductor with a T c of 39 K i 55 ' 56 Based on these 
experiments, we set the small and large superconducting 
gaps as, 



A° =A° =A S 
A 

Q'3 



Q.2 

A 

Q!4 



(+: for s 



±A L , 
-: for s±) 



(25) 



Hereafter, we employ these values as 2Al/T c = 5 and 
Al/As = 3, which are approximately satisfied in vari- 
ous Co-doped 122 systems according to the specific heat 
measurements^ 



Figure [2] (a) shows the temperature dependence of 
7 = ~//V 2 s T 2 on the s ++ -wave superconducting state 
at various frequencies. The T dependence of 7 is small 
in the normal state (T > T c ), but 7 rapidly decreases 
with decreasing T in the superconducting state (T < T c ). 
Figure^] (b) shows the oj dependence of 7. Two vertical 
dotted lines show the small and large superconducting 
gaps at zero temperature As and Al, respectively. In 
the normal state, 7 is large due to the strong correla- 
tion. However, 7 is strongly suppressed in the supercon- 
ducting state since the inelastic damping 7 is reduced 
as the superconducting gap opens. It is satisfied that 
j(u))\t=o = for oj < 3A S = A L . 

Figures [3] (a) and (b) show the T and oj dependence of 
7 on the s±-wave superconducting state, respectively. As 
is the case in s ++ , 7 is strongly suppressed for T < T c . 
The obtained quasi-particle damping for s±-wave state 
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FIG. 2. (a) The temperature dependences of 7 = j/V^ a Tc 
in the s++-wave state with 2Al/T c = 5 and Al/As = 3. 
Solid, dashed, short-dashed and dotted lines corresponds to 
oj — 0, As, Al and 2Al, respectively, (b) uj dependences of 
7 in the s++-wave state for various T. Solid, dashed-dotted, 
dashed and dotted lines represent T = 0, 0.5T C , 0.9T C and T c , 
respectively. Vertical dotted lines show uj = Al and As. 
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FIG. 3. (a) The T dependences of the 7 in the s±-wave state 
with 2Al/T c = 5 and Al/As = 3. (b) The ui dependences of 
the 7 in the s±-wave state. 



is similar to the s ++ -wave state. 



D. Quasiparticle Damping 

In this subsection, we estimate the absolute value of 
the inelastic damping rate 7(0) in the normal state, and 
thus yields the effective interaction V c r , from the exper- 
imentally observed resistivity. From the Nakano-Kubo 
formula, the conductivity is given by, 



g2 ^ 1 dk n ^ fc| 



4ttc-W FSq "27(0)' 



(26) 



TABLE I. T c and other parameters of Ba(Fei_ a: Co :E )2As2 for 
various doping rate 1.— 7(0) is estimated by fitting the ex- 
perimental data and V e ff is obtained by 7(0) and N(0) = 
0.66 eV _1 . We note that we obtain n = 1.0 and a = 0.63 
/ificm/K for x — 0.08 according to the measurement by Sefat 
et aZ.— In this case, we obtain V e g = 19.9 eV. In this table, 
the units of p and a are fiflcm and /if2cm/K n , respectively. 





Ba(Fci_ 


-kCo x ) 2 As2 






X 


0.08 


0.12 


0.14 


0.20 


Tc 


24 K 


17 K 


11 K 




P(0) 


73 


65 


65 


57 


p(To) - p(0) 


11.9 


4.6 


1.3 




a 


0.495 


0.153 


0.045 


0.0011 


n 


1.0 


1.2 


1.4 


2.0 


7(0)|t=t c 


0.0059 


0.0023 


0.00065 






(= 69 K) 


(= 27 K) 


(= 7.5 K) 






17.7 


15.5 


12.7 


4.3 



and 7(0), 27 
p(T) - p(0) 



'0.00207(0) [/^ficm] for c = 6 A 
0.00287(0) [nflcm] for c = 8 A 



(27) 



The interlayer spacing c « 6 A and c ~ 8 A are corre- 
sponding to the 122 and 1111 systems, respectively. On 
the other hand, we solve Eq. (|21l) as a normal state 
(Im/°' R (w) = 0) and the inelastic damping rate 7(0) is 
represented as 



7(0)|r=3 



*-V? s N{0fTl 



(28) 



Then, the effective interaction V e e is derived by compar- 
ing Eq. (g7]) with Eq. 

In Table U we show T c and other parameters of 
Ba(Fei_ a; Co :E )2As2 for various doping rate x. The resis- 
tivity due to inelastic scattering, p{T) — p(0), is estimated 
by fitting the experimental data below 150 K 59 with 



p(T) = p(0)+aT» 



(29) 



where p(0), a and n are free parameters. The obtained 
7 = 69 K for the optimally doped compound {x = 0.08) 
is larger than the transition temperature T c — 24 K. On 
the other hand, 7 decreases rapidly by carrier doping, and 
7 = 7.5 K becomes smaller than T c = 11 K for the over 
doped compound (x = 0.14). The estimated effective 
interaction V e s is 17.7 eV for the optimally doped system, 
and 12.2 eV for the over doped system. V e s is large 
as compared with the bare on-site Coulomb interaction 
U = 2 ~ 3 eV for 3c? electrons on a Fe atorn^ because of 
the spin and charge Stoner enhancement that give large 
spin and orbital fluctuations. 



where 7(0) is the "unrenormalized" damping at zero en- 
ergy, and v a> k is the Fermi velocity at k on the a th Fermi 
surface. Here, we neglect the current vertex correction 
since it is not important for the diagonal conductivity^. 
7(0) is derived from the theoretical relation between p(T) 



III. RESULTS 

In Figs. 0] and we show the numerical results for 
T dependence of the nuclear spin relaxation rate for op- 
timally and over doping states, respectively. According 
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(a) s++ v eH =17.7[eV] 
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(b) s ± V eff =17.7[eV] 
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FIG. 4. The normalized nuclear relaxation rate 
(1/TiT)/(1/TiT)\t=t c for optimally doping state 
(T c = 0.0025 eV and V cfi = 17.7 eV) in the (a) s++ 
and (b) s±-wave states as a function of reduced temperature 
T/T c . Solid and dashed lines represent m*/m = 2 and 
m*/m = 3, respectively. 



to Table HI effective potential V^ff = 17.7 eV for opti- 
mally doped systems and 12.7 eV for over doped sys- 
tems. The mass enhancement factors are m* /m = 2^3 
for optimum doped systems and m*/m = 1 ~ 2 for over 
doped system4^ Figures H (a) and (b) show l/T x T 
in the s ++ and s±-wave states for optimally doping 
with T c = 0.0025 eV and V c g = 17.7 eV. Solid and 
dashed lines represent m*/m — 2 and m*/m = 3, re- 
spectively. In this case, the Hebel-Slichter coherence 
peak is suppressed even in both the s ++ and s±-wave 
states due to the strong inelastic quasi-particle damping 
7*(0)| Tc = 1.0~1.5T c . 

Figures [5] (a) and (b) show \/T\T for over-doped sys- 
tems (T c = 0.0010 eV and V cff = 12.7 eV) in s ++ and s±- 
wave states, respectively. The coherence peak increase 
for smaller V e s since 7 is proportional to V^ s . Since the 
coherence peak is suppressed by 7* = (m/m*)j, it is 
more restored by larger m*/m. Small coherence peak is 
recognized in over-doped s ++ -wave state for m* /m = 2. 
The insets of Figure [S] show l/T x T for T c = 0.0005 eV 



(a) s ++ v e „ =12.7[eV] 




FIG. 5. The normalized nuclear relaxation rate 

(1/TiT)/(1/TiT)| t =t c for over doping state (T c = 0.0010 eV 
and V c f! — 12.7 eV) in the (a) s++ and (b) s±-wave states as 
a function of reduced temperature T/T c . Solid and dashed 
lines represent m*/m = l and m* /m = 2, respectively. 
1/TxT for T c = 0.0005 eV and Kff = 12.7 eV are shown in 
the insets. 



and V e fi — 12.7 eV. Since 7 decreases in proportion to T c 2 , 
small coherence peak appears in both s++ and s±-wave 
states. These results indicate that the Hebel-Slichter co- 
herence peak on \/T\T might be observed for a sample 
with T c ~ 5 K, as reported in Ref. [3^. 

We discuss the Hebel-Slichter peak in both s ++ and 
s ± -wave states. We show the oj dependences of X(w) 
(denned in Eq. ([IB])). iV n (u;) and A^ a (o;) for over doped 
systems (T c = 0.0005 eV, V eS = 12.7 eV and m*/m = 2) 
in Figs. [6J (a), (b) and (c), respectively. The integral 
of X(w) = tt(N d (uj) 2 + iV a (u;) 2 )(-§|j) yields l/T t T. 
Two peaks in X{oj) at T = 0.9T C correspond two su- 
perconducting gaps As and ±Al- These peaks yields 
the Hebel-Slichter peak in the inset of Fig. [5] For s±- 
wave state, the sign of 7V a reserves at oj ~ (Al + As)/2. 
Because of the relation |A fa ++ | ~ |iVJ | in the present 
case Al : As =3:1, similar Hebel-Slichter peak appear 
in both superconducting states. On the other hand, for 
optimum doped systems (T c = 0.0025 eV, V c g = 17.7 eV 



7 



V eB =12.7[eV] 
r Ci =0.0005[eV] 
*=2 

r/r 
1.0 - 




(d) l/ etf =i7.7[eV] 



7" c =0.0025[eV] 



<a=A. 



(e) 



(f) 



m 

S (7=0.97"„) 

®=A L ( T=0 . 9T i 777- 

0.9 (. 
0.9 (. 
1.0 



4 



2 



FIG. 6. (a) The uj dependence of X (uj) / X (0)|t=t c in over 
doped systems for T c = 0.0005 eV, V cB = 12.7 eV and 
m*/m = 2. Solid, dashed and dotted lines represent s++ 
and s±-wave states (T = 0.9T C ) and normal state (T = T c ), 
respectively. (b) A n (w) = -±Im 5 R (o;). (c) 7V a (aj) = 
-±Im/ R (o;). (d)-(f) X(w), 7V n and 7V a in optimum doped 
systems for T c = 0.0025, Kff = 17.7 eV and m*/m = 3. 
Vertical dotted lines show uj = As and Al at T — 0.9T C . 



and m* /m — 3), these peaks in X(ui) are suppressed due 
to strong 7 at T = 0.9T C as shown in Figs. [6](d) - (f). In 
this case, the Hebel-Slichter peak disappears as shown in 
Fig. 1 



IV. SUMMARY AND DISCUSSION 

To summarize, we investigate the nuclear magnetic re- 
laxation rate 1/T\T, especially paying attention to the 
presence of the coherence peak for both s++ and s±- 
wave states base on the five orbital model. The inelastic 
quasi-particle damping rate 7 at T c is estimated from 
the experimental results of the resistivity, and T- and in- 
dependences of 7 are calculated using the second order 
perturbation theory. 

Using parameters for optimally doped systems, the 
Hebel-Slichter coherence peak on \fT{F is suppressed 
due to the strong 7*(0)|t c ~ T c even in the s ++ -wave 
state for T > 10 K. This result is consistent with pre- 
vious strong coupling theoriesi 36 i 37 On the other hand, 
the relation 7*(0)|t c <C T c is expected in heavily over 
doped systems with T c ~ 5 K: In this case, coherence 
peak may appear for both pairing states. Note that the 
tiny coherence peak in Fig. [5] (b) in the s-t-wave state 
grows comparable to that in [5] (a) in the s++-wave state 
by halving the value of 7*(0). Thus, the condition for the 
appearance of coherence peak is similar for both s ++ and 
s±-wave states, so it is difficult to discriminate between 
these pairing states from the present NMR experimental 
data. 

In Refs. and 0, the authors measured \jT\ 



and discussed the effect of 7 on the coherence 
peak. They assumed the functional form 7(0;) = 
7(0)|t=t c expL4 (T/T c — 1)] (uj dependence was ne- 
glected), and chose parameters A — 5 and 7(0) \t c = 3T C , 
and found that the coherence peak disappears even in 
the s ++ -wave state. The large damping 7(0) |r c = 3T C is 
comparable to our estimation (1.0 ~ 1.5T C ) for optimum 
doping systems, so their analysis is consistent with the 
present study. 

In this paper, we discussed only Co-doped Bal22 sys- 
tems, since reliable resistivity data in single crystals 
are available. The coherence peak is also absent in 
LiFeAs (c = 6.36 A, T c = 17 K) ?2 . In this com- 
pound, the resistivity is fitted as p(T) = po + aT 2 and 
a = 0.022 /iOcm/K 2 in Ref. H Then, the obtained 
effective potential is V e g — 15.9 meV, which is com- 
parable to the value of optimally doped Ba(Fe,Co)2As2 
(17.5 meV). Therefore, if we assume the gap structure 
given by ARPES measurements (A L - A s — 3 meV)£^£ 
the coherence peak disappears even for the s++-wave 
state. However, since the size of each hole-pocket in 
LiFeAsH££ are different from that in Ba(Fe,Co) 2 As 2 , 
more realistic tight-binding model for LiFeAs would be 
required for a quantitative analysis. 

In LaFeAsOi-^Fj,, we cannot estimate 7 quantita- 
tively because of the lack of resistivity data in single crys- 
tals. For a qualitative analysis, however, we can roughly 
estimated the single crystal resistivity of LaFcAsOi-^Fa; 
from the poly crystal resistivity in Ref. l66l by multiply- 
ing the factor 1/3, as discussed in Ref. [Ill Using this 
method, we can fit them by using p(T) = p(0) + aT 2 
and obtain 7(0) | Tc ~ 0.0023 eV and V cff - 10.2 eV for 
LaFeAsO0.89F0.11 (T c — 26 K). Using these parameters, 
the Hebel-Slichter coherence peak for optimally doped 
Lallll systems are also suppressed in both s++ and s±- 
wave states. 

In both LiFeAs and Lallll, conventional Fermi liq- 
uid resistivity (p oc T 2 ) is observed. This fact would 
be consistent with the small spin fluctuations in these 
compounds. However, 7 at T c is large enough to sup- 
press the coherence peak. Thus, both compounds are 
strongly correlated Fermi liquids due to relatively local 
fluctuations. In Smllll and Ndllll (T c > 50 K), the 
T-linear-type behavior of resistivity is observed, indicat- 
ing the enhancement of fluctuations as AS4 tetrahedron 
approaches to a regular tetrahedron^ 

We stress that the inelastic scattering 7 also plays 
important roles in the neutron scattering spectrum: 
In neutron magnetic scattering measurements, broad 
" resonance- like" peak structure is observed in many FeAs 
superconductors below T c , and this fact had been fre- 
quently ascribed to the resonance due to the gap sign 
change. However, the "resonance condition" is not surely 
confirmed since it is difficult to determine the gap size 
accurately. The resonance condition in FeAs supercon- 
ductor is uj rQS < Al + Ag, where w rcs denotes the peak 
energy of neutron scattering. In the following, we write 
down the experimental values of w ros , Al and As- 
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(A) BaFei.85Coo.i5As2: uj ics = 10 meV by neutron^ 

* Specific heati^i 

Al = 5 meV, A s = 2 meV: A L + A s = 7 meV 

* Penetration depths 

A L = 6.1 meV, A s = 2.3 meV: A L + A s = 8.4 meV 

* ARPES)^ 

A L = 6.6 meV, A s = 5 meV: A L + A s = 11.6 meV 

(B) FeTeo.6Se . 4 : w res = 7 meV by neutron^ 

* ARPESiii 

A L = 4.2 meV, A s = 2.5 meV: A L + A s = 6.7 meV. 

In case (A), the resonance condition oj Tes < Al + As 
is satisfied only by ARPES data. In case (B), although 
ARPES tends to report larger gap size in FeAs supercon- 
ductors, the resonance condition is still unclear. 

Theoretically, we have revealed that broad hump struc- 
ture at uj ~ Al + As can appear even in the s++-wave 
state, not due to the resonance but due to strong sup- 
pression of inelastic quasi-particle scattering rate 7(0;) 
for < 3 A (" dissipation- less mechanism" 26 i 27 ) . Thus, 
both the absence of the coherence peak in 1/T±T as well 



as the hump structure in the neutron scattering spectrum 
are explained by the same many-body effect — inelastic 
quasi-particle scattering — using similar inelastic scat- 
tering rate 7(0) evaluated from the resistivity. 

Finally, we comment that large inelastic scattering 7 at 
T c in Fe based superconductors also works as the depair- 
ing effect. In the random-phase-approximation (RPA), 
in which the depairing effect due to inelastic scattering is 
neglected, the obtained T c is about 200 K£ii£ However, 
in the fluctuation-exchange (FLEX) approximation, the 
obtained T c is strongly suppressed (below 50 K) due to 
the the depairing effect (self-energy correction)^ 
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